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Abstract 

The general linear quaternion function of degree one is a sum of terms with quaternion coefficients 
on the left and right. The paper considers the canonic form of such a function, and builds on the 
recent work of Todd Ell, who has shown that any such function may be represented using at most 
four quaternion coefficients. In this paper, a new and simple method is presented for obtaining these 
coefficients numerically using a matrix approach which also gives an alternative proof of the canonic 
forms. 



1 Introduction 

In a recent paper [ ] , Todd Ell has shown that a linear quaternion function of the form: 

p 

fil) ^^^pI'T^P (1) 

where the coefficients mp and n^, and the variable q, are quaternions, can be expressed in a canonic 
form using a maximum of four quaternion coefficients, irrespective of the number of terms P. Ell's 
paper gives a method to obtain the four coefficients of the canonic form algebraically without explicit 
decomposition of the original quaternion coefficients nip and Hp into their four components. This paper 
does not improve on that result, but it does present a simple method for obtaining the four coefficients 
of the canonic form numerically. The canonic form is^: 

f{q) = Aq + Bqi + Cqj + Dqk (2) 

Subsequent to the publication of this result by Todd Ell, a paper by Littlewood in 1931 was found to 
contain this canonic form [.''] although Littlewood gave no justification for his statement that it was 
'general' (assumed to mean 'canonic') nor a method to obtain the four coefficients from the general case 
as in equation 1. No earlier papers are currently known that consider this form. In a later paper with 
Richardson, the following equation appeared [ ; , page 333]: 

n n 
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The context of this equation is non-commutative algebras in general. In the specific case of quaternions 
Co = 1, ei = i, 62 — j and 63 — fc, and the result given in equation 3 in the next section is a special case 
of Littlewood and Richardson's expression. 
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^It is also possible to obtain a similar form with the quaternion coefficients on the right and the quaternion operators 
i, j and fc on the left, as was shown in [ ]. 



2 Preliminaries 



Consider a single term in equation 1: mpqup where we write = Wp + Xpi + ypj + Zpk and Up = 
w'p + x'pi + y'pj + z'pk. The product may be expanded as follows: 

rupq Up = {wp + Xpi + ypj + Zpk)q{w'p + x'pi + y'pj + z'pk) 

Wpw'p q + WpX'p qi + Wpy'p qj + Wpz'p qk + 
^ Xpw'piq + Xpx'piqi + Xpy'piqj + Xpz'^iqk + (3) 
~ Vpw'pjq + ypx'pjqi + ypy'pjqj + Vpz'pjqk + 

Zpw'pkq + Zpx'pkqi + Zpy'pkqj + Zpz'pkqk + 

We can immediately see that the full function in equation 1 can be written as a sum of sixteen summations 
over the products of components of the left and right quaternion coefficients: 



p 




.f{'l) = Y^p1"'P^ /p^ \ \ /p^ \ /p^ \ 



This expansion shows immediately that at least two canonic forms exist, as already shown in [ ]. Group- 
ing the four summations within the same column, each column sums to give a quaternion coefficient on 
the left of q: 

f{q) = Aq + Bqi + Cqj + Dqk 

An alternative is to group the four summations within the same row, where each row sums to give a 
quaternion coefficient on the right of q: 

f{q) = qA' + iqB' + jqC + kqD' 

Now, a significant question is whether this is all, or whether there are additional canonic forms derivable 
from this expansion. More accurately, are there additional canonic forms with 4 quaternion coefficients 
(since a canonic form with more than 4 coefficients would be sub-optimal)? The next section recasts the 
problem in matrix form in order to show how this question might be answered. 

In case this is thought by readers to be a trivial issue, Ljudmila Meister [5, §3.5] proposed the following 
as a general linear form: f{q) — Aq + qB + CqD. This has the correct number of quaternion coefficients 
(four, the same as the canonic form given by Littlewood in 1931 [']) and it looks very plausible as a 
canonic form, since there can be only two terms with a single quaternion coefficient, one on the left 
and one on the right (any others could be trivially combined with these since Aq + Eq = {A + E)q), 
and one double-sided term with two coefficients. In fact, using the matrix formulation presented in the 
next section, it is argued that this form cannot be canonic. It is very hard to show this by an algebraic 
argument, but the matrix formulation presented in the next section makes it relatively simple. 



3 Matrix and outer product formulation 



The coefficients in equation 3 can be expressed as the following matrix, which is the outer product [ ] of 
two vectors representing the left and right quaternion coefficients and Up. It follows from the fact 
that this matrix is the outer product of two vectors that it must be of rank 1 [' ] . 
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Similarly, we can express the coefficients in equation 4 as the following matrix. It follows from the fact 
that this matrix is 4 x 4 that it must be of rank 4 or less. 
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M may be decomposed into the sum of four matrices of rank 1, each of which can be factored into an 
outer product of two vectors as in equation 5. It follows that, for any number of terms P, equation 1 
can be expressed as at most the sum of four terms, that is P = 4. The four rank 1 matrices which are 
summed to make M can be obtained from a singular value decomposition [ ], which can also yield the 
vectors whose outer product gives one of the components of M. (Products of other pairs of columns will 
be zero because of orthogonality.) This can be done as follows: 



where U and V are orthogonal, and Z is diagonal. Corresponding columns of matrices L and V are vectors 
whose outer products give the four matrices which sum to M. These vectors are not unique, although 
the singular values in Z are. 

Clearly, a general linear quaternion function of the form in equation 1 with an arbitrary number of 
terms greater than 4 will yield a matrix M of rank 4. It follows therefore that any such function can 
be expressed using at most 4 terms, where the left and right quaternion coefficients in each term are 
obtained from the outer product factorisation of rank 1 matrices composing M (these rank 1 matrices 
can be obtained from a singular- value decomposition of M): 

/((?) = AqE + BqF + CqG + DqH 

However, this is clearly not optimal, since it has eight quaternion coefficients, double the number needed 
in the canonic form given in equation 2. However, by decomposing matrix M into the sum of four matrices 
each containing one column of M, like this: 



(7) 



we obtain the canonic form f{q) = Aq + Bqi + Cqj + Dqk in which the right-side quaternion coefficients 
reduce to the degenerate set 1, i,j,k and the left-side coefficients are given hy A = mn + m2ii + m^ij + 
muk and so on. (We could do the same with the rows of M, in order to obtain the canonic form with 
1, i, j and k on the left.) 

We now return to the 'general' form proposed by Ljudmila Meister [■">, §3.5] and mentioned in the 
last paragraph of §2. This form is: f{q) = Aq + qB + CqD. It should be clear from the foregoing that 
the three terms in this function would yield a matrix M of rank 3, and since the general case requires 
a matrix of rank 4, it is not canonic. It might appear that adding a fourth term (EqF) would fix the 
problem, since this would increase the rank to 4, but in fact this is not sufficient to yield a canonic 
form. This can be seen as follows. A corresponds to a column of the matrix M, and P to a row. This 
leaves three other columns and rows, which in the general case, will be of rank 3. It follows that two 
double-sided terms fike CqD and EqF are not sufficient to represent the rest of the terms of a general 
linear function, since they cannot yield a component of M of rank greater than 2, and a rank of 3 is 
needed. We are thus led to conclude that a function of the form f{q) — Aq + qB + CqD + EqF + GqH 
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is required, which is again not optimal, since it has 8 quaternion coefficients. In fact, it is not necessary 
for all the terms in this function to be full quaternions, but as we now show, this docs not improve on 
the form in equation 2. 

If we take the first column and the first row for A and B, we would be left with a 3 x 3 matrix, which 
in the general case would be of rank 3. Factorising it into outer products we would need three terms 
with a pure quaternion on each side. Noting also that it is not necessary for both of A and B to be full 
quaternions, a 'canonic' form appears to be: f{q) = Aq + qb + Viqv2 + v^qv^ + v^qv^, where b and the Vi 
are vectors (pure quaternions). The number of real coefficients in this function is 25, 9 more than in the 
canonic forms. However, when factorising the 3x3 matrix within M we can factorise by columns, thus 
reducing the form to f{q) — Aq + qb + viqi + v^qj + v^qk which has 16 real coefficients and is therefore 
canonic, although it is just a variation of the form derived recently by Ell, and published by Littlewood 
in 1931. 

It appears unlikely that there are any other decomposition(s) of M that would yield a canonic form 
with only 16 real coefficients. 

4 Conclusion 

In general, an arbitrary quaternion function of the form given in equation 1 cannot be reduced to less 
than four double-sided terms. That is, the function can always be reduced to a sum of 4 double-sided 
terms, with 8 quaternion coefficients. There are useful special cases where four of these coefficients are 
degenerate either on the left, or on the right (but not mixed) as already shown in [1]. 

The matrix formulation given in this paper provides a simple and systematic way to study the problem 
further, although it appears unlikely that the canonic form in equation 2 can be improved on. 
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